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SUMMARY 


A general procedure for computing the region of influence of a maneuvering 
vehicle is described. Basic differential geometric relations, including the 
use of a general trajectory parameter and the introduction of auxiliary vari- 
ables in the envelope theory are presented. To illustrate the application of 
the method, the destruct region for a maneuvering fighter firing missiles is 
computed. 


INTRODUCTION 

In the course of the flight of an airplane, a missile, or even a space 
vehicle, certain regions in the neighborhood of the flight path are of inter- 
est, usually because of the potential influence that the vehicle has on the 
space in its vicinity. This influence may be of a passive nature, such as 
noise, engine exhaust, or trailing vortices. It may, on the other hand, be 
deliberate in character. The vehicle may, for example, emit missiles or other 
projectiles, radar seeking or jciraming signals, and laser guidance signals. 

In such a situation it is useful to know the extent of the region of 
influence as a function of the vehicle's path and its speed. The shape of 
the region is generally known for constant conditions, that is, for the vehi- 
cle at rest or in straight steady flight. The total region of influence for a 
general flight trajectory is the interior of the envelope of all the instanta- 
neous regions. It would be possible to calculate each region separately and 
then determine the envelope numerically by interpolation. However, it is much 
more efficient to calculate the envelope directly by determining the character- 
istic lines (ref. 1) as described in the present paper. Another advantage of 
this method is that the analytical nature of the envelope equations provides a 
qualitative insight into the effects of various flight parameters on the size 
and shape of the region of influence. 

In the present paper the general procedure is described. The basic equa- 
tions are given, together with scxne specialized analytical results that are 
useful in the application of the method. Two illustrative examples, both 
dealing with fighter maneuvers, are given in some detail. These examples 
treat the intercept region as a function of the trajectory of the attacking 
airplane only. 

The basic technique has already been applied specifically to the sonic 
bom problm (ref. 2) , but the generality of the method and the wide range of 
its potential applications in aerospace problems have apparently not been 
recognized. Some of these applications include the calculation of a missile 
intercept region, a jet noise constant-decibel region, and a region limited by 
the lens resolution of a camera carried by a spaceship. 


SYMBOLS 


a speed of sound, m/sec 

b unit vector perpendicular to both T and n, T x n 

c constant of proportionality (used in eq. (8)), 1/m 

0*1 constant relating missile range increment to vehicle velocity 

(used in eq. (25)), 1/sec 

d total maximum range, m 

d^ basic range for projectile fired from rest, m 

i,j,k unit normal base vectors fixed in space 

I local torsion of flight path, 1/m 

M Mach number 

n function shape parameter (used in eq. (24)) 

n unit normal vector (used in eq. (3)) 

p general trajectory parameter 

R radial variable in spherical coordinate system, m 

r position vector of a point on the periphery of the region of 

influence of a maneuvering vehicle, m 

r^ position vector of vehicle, m 

s arc distance along trajectory, m 

t time, sec 

T unit vector locally tangent to flight path at vehicle location 

V velocity, m/sec 

Vfji total velocity of projectile, m/sec 

u,v general variables 

x,y,z components of r 

0 angle coordinate relative to T vector, deg 

0* radar semicone angle, deg 
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K local trajectory curvature, 1/m 

P local trajectory radius of curvature, 1 /k, m 

X time at which vehicle is at point r^, sec 

(j) angle coordinate, deg 

Superscripts: 

’ derivative with respect to s 

• derivative with respect to time 

Subscripts: 

o conditions when a projectile is fired from rest 

p,u,v,0,T derivative with respect to the variable 

ANALYSIS AND EXAMPLES 
Review of Basic Equations 

In this section, mathonatical expressions are given for the moving trihe- 
dral coordinate system and the trajectory curvature and torsion. The position 
vector of the vehicle at time x is 

rt(X) = xi + yj + zk (1) 

The moving trihedral base vectors are determined as follows: 

If s denotes the length along the trajectory path, the vector T is 
given by 

_ dx _ dy _ dz « 

T = — i + — j + — k (2) 

ds ds ds 


which represents a unit vector tangent to the flight path. Therefore, it is 
perpendicular to its derivative, which can be written in the form 


1 dT 

n = (3) 

K ds 
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where n is a unit vector, and ic is the magnitude of dT/ds and is termed 
the trajectory curvature (ref. 1, p. 13). A third unit vector perpendicular 
to both T and n is defined by 


b = T X n 


( 4 ) 


The three mutually perpendicular unit vectors T, n,_and b define the 
moving trihedral coordinate system (fig. 1). The vector T points in the 
direction of flight, and n points in the direction of the local center 
of curvature of the trajectory. Other quantities which are used later are 
the derivatives of n and b which are given by 


dn « - 

— = -KT + Zb (5) 

ds 


db 

— = -zn (6) 

ds 


In equations (5) and (6), which together with equation (3) are known 
as the Frenet-Serret formulas, z is the local trajectory torsion (ref. 1, 
pp. 15-18). 

The trajectory parameters s and x are related through the velocity 
which is assumed to be known either as a function of s or x; 


ds 

— = V (7) 

dx 


Example: Fighter Firing Cannon Shells 

To illustrate the application of these equations within a physical con- 
text, consider a fighter that is firing projectiles from a cannon while maneu- 
vering. The projectiles are fired in the (local) direction of flight, that is, 
in the direction of T. Their deceleration due to aerodynamic drag is pro- 
portional to the square of the speed: 


Vrp = -CVip^ 


(8) 
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Thus, 


^ Vqi dViji 

'^V+Vo 



- c dt 


vrtiich gives for the projectile velocity at time t. 


V + Vo 

Vt = 

1 + c(V + Vq) ( t - T) 


(8a) 


and for the distance traveled from the point of firing. 


d = - In [l + c(V + Vq) ( t - t) ] 
c 


(9) 


If the projectiles are fired over the period t«| < T < t, then the dis- 
tribution of the projectiles in space at time t is given by the vector 
expression 


r = rt(T) + d(t,T) T(T) (tT < T < t) (10) 

where d is given by equation (9) . Figure 2 shows an example of such a dis- 
tribution for an airplane performing a turning maneuver. (In the example, 
c = 0.0676/m, V = 183 m/sec, Vq = 183 m/sec, and p = 1219 m.) 

If it is assumed that the projectile loses its effectiveness after its 
velocity decreases to a certain value, then by inserting this value in the left 
side of equation (8a) and eliminating (t - i) from equations (8a) and (9) , an 
effective range can be calculated. With this range, an alternate problem would 
be to compute the entire "destruct" surface swept out by the shells fired over 
a period of time, allowing for each shell its full range as limited by effective 
velocity. For this problem, the surface is bounded on one edge by the flight 
trajectory rt('^) and on the other by the line 


r = rt('t) + d(V)T (11 ) 


An example of this type of problem is shown in figure 3 for the same type of 
maneuver used for the example of figure 2. In this example, the projectile 
effective range is taken to be d =* 1830 m. 
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In these examples, simplified mathematical models have been assumed 
for purposes of illustration. They could be made more sophisticated with 
some increase in mathematical complexity. For example, the velocity equa- 
tion (8a) could be written as a vector equation with a component in the 
ic-direction resulting from gravitational attraction, and components taken 
in the T,n plane because, in performing a turning maneuver, the direction 
of firing is not exactly in the T-direction. 


Form of Trajectory Equations 

The preceding examples require the differentiation of the flight path 
vector r^ in order to compute the unit tangent vector T (eq. (2)). 

More complex problems require higher derivatives of r^. It is generally 
advantageous to express the flight trajectory in analytic form so that the 
derivatives can be calculated analytically. Analytic expressions would nor- 
mally be used in any hypothetical situation, but even if the trajectory repre- 
sented an actual flight, with locations determined at specific times, it would 
be preferable to fit the trajectory with an analytic curve. 

However, it is not a simple matter to express most curves in terms of a 
specific parameter, such as the arc-length parameter s. Circular and elliptic 
arcs, for example, are more naturally expressed in terms of an angle variable. 
Thus, if the trajectory coordinates are expressed in terms of a general param- 
eter p, then the trajectory vector is of the form 


?t(P) = x(p)I + y(p) j + z(p)k 


and 


T = p' (Xpi + ypj + Zpk) 


( 12 ) 


v^ere 


P' 


1 



(13) 


Also 


dT ^ 

— = P"rtp + P'^rtpp 
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which, when written in component form, is 


dT - _ 

— = (P"xp + p'2xpp)i + 


(P"Yp + P'^pp) j + (P"Zp 


+ P'2zpp)k 


(14) 


where 


P" = P* 


dp- 

dp 


or 


„ ^ _ XpXpp + YpYpp + ZpZpp 

(Xp2 + yp2 + Zp2)2 


Thus, performing all calculations in terms of the_parameter p, one can 
calculate T from equations (12) and (13), and n and K from equa- 
tions (3), (14), (13), and (15). 


Envelope Theory 

The region of interest that is associated with a maneuvering vehicle is 
usually not a simple set of vector lines, as in the example of the cannon 
shells fired from a fighter, but rather is a three-dimensional region. Sup- 
pose, for example, the fighters were firing missiles rather than cannon shells. 
Since the missiles can be fired over a certain cone angle range, there is, 
associated with each point of the trajectory, a destruct region within which 
a target can be destroyed. Then, to determine the boundary of the destruct 
region associated with the entire trajectory, one must determine the envelope 
of all the surfaces that bound these individual regions. Somewhat similar 
problems involving the envelope of a family of surfaces arise in calculations 
of subsonic and supersonic airplane noise. 

The theory of calculating the envelope of a one-parameter family of sur- 
faces is given in reference 1 (pp. 162-168) . If the parameter is p and the 
family of surfaces is given in the form 


F(x,y,z,p) = 0 


(16) 
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then, when the equation 


Pp(x,y,z,p) = 0 


( 17 ) 


is solved simultaneously with equation (16) at a given value of p, the result- 
ing solution determines a characteristic line. The set of these characteristic 
lines, determined as p varies, defines the envelope of the surfaces described 
by equation (16). 

A simple example is the family of sj^eres representing the periphery of 
the pressure disturbance generated by a moving point source (see ref. 2) : 


P(x,y,z,T) = (r - rt) * (r - r^-) - a2(t - T) 2 = q (18) 

Differentiating with respect to the parameter t yields 
-2VT • (r - rt) + 2a2(t - T) = 0 


or. 


_ - - a(t - T) 

(r - rt) • T * 


This is the equation of a plane which, together with the sphere equation (18), 
determines a circle (at supersonic speeds) for any given value of t. At 
time t, the set of these circles formed as t varies from 0 to t comprise 
the Mach surface of the moving source. 

Por this illustration, the basic region, being spherical, is readily 
described in the form of equation (16). But other problems, as for example 
those involving symmetry about just one axis, involve surfaces that are not 
easily expressed analytically in this form. A more versatile expression for 
the equation of the basic surface is 


r * r(u,v,p) (19) 


where u and v represent auxiliary variables or new coordinates. With this 
formulation the characteristic line is obtained by solving equation (1 9) simul- 
taneously with the ec[uation 


■ r5.. 
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( 20 ) 


(See ref. 1 , p. 168.) 

A common case of interest in trajectory problems is that of a region that 
is symmetric about the direction of motion , that is, about the T axis. The 
geometry for this situation is depicted in figure 4. According to the figure 


r =s r^ + R cos 0 T + R sin 9 cos (j) n + R sin 0 sin b (21) 

where, for the axisymmetric region, R is independent of (|); 

R = R(0,T) (21a) 

Thus, if the trajectory parameter is taken to be T, 
r = r(0,(|),T) 


which, if (0,(|),T) are substituted respectively for (u,v,p) , is in the form of 
equation (19). Thus equation (20) yields 


R0 COS 6 - R sin 0 (Rq sin 0 + R cos 6) cos (j) (R© sin 0 + R cos 0) sin (j) 

0 -R sin 9 sin (J) R sin 0 cos (j) 

I V(1 - KR sin 0 cos 4>) + Rj cos 0 VR(< cos 0 - i sin 0 sin <J)) + R^- sin 0 cos (VRz cos ({> + R^- sin (})) sin 0 


> = 0 


Expanding this determinant and collecting terms yields the equation 


rRt^ 

K7 sin 9 [-RR 9 ic cos (|) + R 0 sin 9 + R cos 0 + 1 = 0 


( 22 ) 


It is interesting that, in evaluating the determinant, the torsion terms sub- 
tract out. 

For a straight trajectory k = 0, and the first term in equation (22) 
drops out. The resulting equation is independent of cj), and consequently the 
characteristic line is a circle. When the trajectory is curved, equation (22) 
can be written in the form 
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(23) 


R sin 6 cos 4> = P sin 


9 (sin 0 


R „ ' R 

+ — cos 0 + 1 

R0 R0 V y 


unless the range of 0 includes values for which R0 - 0. The left side of 
equation (23) is the component of the characteristic vector in the direction of 
n (eq. (21)). The corresponding ccanponent in the T-direction is R(9) cos 0. 
These two components and the vector length R, obtained frc»n equation (21a), 
are sufficient to determine the component in the b-direction. Thus, in this 
case, the characteristic line can be plotted as a function of 9 without com- 
puting the corresponding ({> values explicitly. 

Example: Missile Destruct Region 

An illustration of the use of these equations is provided by the previ- 
ously mentioned problen of determining the region within which a target can be 
destroyed by air-to-air missiles fired by a maneuvering airplane. For this 
example the basic region is considered to be determined by the radar visibility 
region and the missile range. 

The radar visibility region is a function of the radar sweep azimuth and 
elevation angles and may be asymmetric and time-varying in form. However, for 
simplicity in the illustration, it is approximated by a fixed-angle cone. The 
missile range is determined as a basic range cIq, when fired frcxn rest, together 
with a component Ad imparted to it by the forward motion of the airplane. The 
total range is therefore a maximum in the direction of flight and falls off 
scxnewhat if the missile is fired at an angle to this direction. 

The total destruct region, then, is limited laterally by the radar scan 
and longitudinally by the missile range, as depicted diagrammatically in fig- 
ure 5 (a) . This region will be approximated by the region represented by the 
analytic expression 


R = d(T) cos n0 = [dQ + Ad(T)] cos n0 


(24) 


where d is the maximum range and the parameter n is chosen so that the ray 
to the point of maximum diameter matches the radar cone angle. (See fig. 5(b).) 
The point of maximum diameter is determined by differentiating the expression 
for the diameter 


R sin 9 = d cos n0 sin 9 


and setting the derivative equal to zero: 


d(cos 9 cos n9 - n sin n9 sin 9) = 0 
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Thus, if 0* is the radar semicone angle, the value of n satisfying 


n = cot 9* cot n0* 


is to be inserted in the approximating equation (24) . If the simple linear 
estimate 


Ad(T) * ctV 


( 25 ) 


is assumed, then equations (23) and (24) yield 


R sin 0 cos (t> = P sin 9 


cot n9/' CiVt cos n0 

sin 0 - cos 9 + 

n V V , 


(26) 


The remaining components of the characteristic vector can be obtained as 
explained following equation (23) . 

Figure 6 shows, in perspective view, an example for a radar angle of 16® 
and a maximum range of 3.11 km. The vehicle trajectory has an increasing 
curvature, giving rise to an increasing elongation of the characteristic 
curves and, consequently, of the envelope. 

Thus the intercept region for a missile in straight flight is confined to 
a relatively narrow cylinder, but the lateral extent of the intercept region 
for an airplane on a helical trajectory is significantly greater. 

The shape of the characteristic lines is strongly sensitive to the trajec- 
tory curvature when the basic region given by equation (24) is highly elongated. 
This dependence can be observed qualitatively by studying the envelopes for 
various values of the parameter n in equation (24) . For n = 1 , the region 
is a sphere and the characteristic lines are circles. 

Figure 7 shows, for a slightly curved trajectory (P = 3.05 km), the enve- 
lope that results for values of n = 1.5, 2.0, and 3.0 corresponding, respec- 
tively, to radar angles of 31.5®, 24.1®, and 16.3®. It is seen that, as the 
radar angle 9* decreases, the characteristic lines and the envelope become 
highly elongated. 


CONCLUDING REMARKS 

A general procedure for computing a region of influence associated with a 
maneuvering vehicle has been described. The method described for calculating 
the envelope of the individual constant- state regions directly from the ana- 
lytic expression for the characteristic lines has significant advantages over 
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the strictly numerical approach. It is more efficient computationallyr and 
the analytic form of the characteristic equation provides a qualitative insight 
into the effect of the various flight parameters on the size and shape of the 
region of influence. Application of the theory was illustrated by computation 
of the destruct region for a fighter firing cannon shells, as an example of 
a simple vector line surface, and the destruct region for a fighter firing 
missiles, as an example of the application of the envelope theory. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 20, 1980 
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Figure 1 Fixed and moving coordinate systems. 


Projectile path 



Diagrammatic view 



Computer plot 

Figure 2.- Example of computing distribution at fixed time 
of projectiles fired from turning airplane. 





Computer plot 

Figure 3.- Destruct surface swept out by projectiles fired fran 
turning airplane (perspective view) . 
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Boundary determined 



(a) Potential missile effectiveness region as determined by 
radar sweep angle and missile range. 



(b) Analytic approximation to missile effectiveness 
region in the form d cos n0. 

Figure 5.- Physical and analytic approximations to 
missile effectiveness region. 



Flight path 
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Airplane at time 

Airplane at time T2 

Missile effectiveness region 
for airplane at 

MissUe effectiveness region 
for airplane at T2 


Characteristic line of 
envelope, tangent to both 
regions 



Set of characteristic 
lines constituting 
envelope 


(a) Concept schematic. (b) Computed example. 

Figure 6.- Section of missile effectiveness envelope for airplane 
performing turn of increasing curvature. 



(a) 0* = 31.5°. 


(b) 9* = 24.1°. 


(c) 9* = 16.3°. 


Figure 7.- Comparison of envelope shapes for increasing elongation of basic effectiveness 
region (decreasing radar cone angle) , for slightly curved flight path. 
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